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ABSTRACT 

The origin of the magnetic fields observed in some intermediate mass and high mass main 
sequence stars is still a matter of vigorous debate. The favoured hypothesis is a fossil field 
origin, in which the observed fields are the condensed remnants of magnetic fields present 
in the original molecular cloud from which the stars formed. According to this theory a few 
percent of the PMS Herbig Ae/Be star should be magnetic with a magnetic topology similar 
to that of main sequence intermediate-mass stars. 

After our recent discovery of four magnetic Herbig stars, we have decided to study in 
detail one of them, HD 200775, to determine if its magnetic topology is similar to that of the 
main sequence magnetic stars. With this aim, we monitored this star in Stokes / and V over 
more than two years, using the new spectropolarimeters ESPaDOnS at CFHT, and Narval at 
TBL. 

By analysing the intensity spectrum we find that HD 200775 is a double-lined spectro- 
scopic binary system, whose secondary seems similar, in temperature, to the primary. We have 
carefully compared the observed spectrum to a synthetic one, and we found no evidence of 
abundance anomalies in its spectrum. We infer the luminosity ratio of the components from 
the Stokes / profiles. Then, using the temperature and luminosity of HD 200775 found in the 
literature, we estimate the age, the mass and the radius of both components from their HR 
diagram positions. From our measurements of the radial velocities of both stars we determine 
the ephemeris and the orbital parameters of the system. 

A Stokes V Zeeman signature is clearly visible in most of the Least Square Deconvolution 
profiles and varies on a timescale on the order of one day. We have fitted the 30 profiles 
simultaneously, using a x 2 minimisation method, with a centered and a decentered-dipole 
model. The best-fit model is obtained with a reduced^ 2 = 1 .0 and provides a rotation period of 
4.3281+0.0010d, an inclination angle of 60+11°, and a magnetic obliquity angle/? = 125+8°. 
The polar strength of the magnetic dipole field is 1000 + 150 G, which is decentered by 
0.05 + 0.04 R„ from the center of the star. The derived magnetic field model is qualitatively 
identical to those commonly observed in the Ap/Bp stars. 

Our determination of the inclination of the rotation axis leads to a radius of the primary 
which is smaller than that derived from the HR diagram position. This can be explained by a 
larger intrinsic luminosity of the secondary relative to the primary, due to a larger circumstellar 
extinction of the secondary relative to the primary. 

Key words: Stars : magnetic field - Stars: pre-main-sequence - Stars: binaries: spectroscopic 
- Star: individual: HD 200775 - Instrumentation : spectropolarimetry. 



Based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council of Canada, 
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Table 1. Log of the observations. Columns 1 and 2 give the UT date and the Heliocentric Julian Date of the observations. Column 3 gives the total exposure 
time. Column 4 gives the peak signal to noise ratio (at ~ 730 nm per ccd pixel) in the spectra and column 5 gives the signal to noise ratio in the LSD Stokes 
V profiles. Column 6 gives the longitudinal magnetic field, column 7 gives the rotation phase derived in Section 6.2, and columns 8 and 9 give the radial 
velocities of both components of the system. Column 10 gives the instrument used. 
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1 INTRODUCTION 

Some main sequence A, B and O stars host strong (~ kG) organised 
magnetic fields. The origin of the magnetic fields of these interme- 
diate and high mass stars is still a matter of debate. The favoured 
theory is the fossil field hypothesis. This theory assumes that the 
magnetic fields observed in these main sequence stars are relics 
of the magnetic fields which existed in the molecular clouds from 
which the stars formed. This theory implies that the remnant inter- 
stellar magnetic field should subsist throughout all the processes of 
formation encountered by the star, from the gravitational collapse 
to the main sequence phase, without being regenerated. According 
to this theory, some pre-main sequence (PMS) stars of intermediate 
mass should host magnetic fields. 

For a long time, this theory conflicted with the general belief 
that all stars pass through a completely convective phase during the 
PMS phase. The turbulent diffusion produced by the convection 
wou ld dissipate the m agnetic field during the Hayashi phase. How- 
ever IPalla & StahleiUl993h calculated the birthline (the locus in the 
HR diagram where stars become optically visible and start the PMS 
phase), showing that the Hayashi phase is considerably reduced for 
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stars between 1.5 M Q and 2M Q and disappears completely for stars 
above 2M Q . A fundamental conclusion of this work was that the 
magnetic field of the intermediate PMS stars can potentially sur- 
vive the PMS phase. 

According to the fossil field hypothesis, some Herbig Ae/Be 
stars (PMS stars of intermediate mass) should be magnetic. Many 
authors have tried to detect a magnetic field in these stars. 
ICatala et al. I dl993l) observed AB Aur using a Zeeman polarime- 
ter placed before the Coude spectrograph of the CFHT They found 
no circular polarisation signal in the Fell 5018 A line a nd they ob- 
tained an upper limit around 1 kG. ICatalaetal.i l ! 1999b tried again 
to detect a magnetic field in AB Aur using the MUSICOS instru- 
ment temporarily installed at the CFHT. They observed no Zee- 
man signature and obtained a lower detection limit around 300 G. 
iDonati et al.l dl997t) observed a large sample of cool and hot stars, 
including two Herbig stars, using the UCLES spectrograph of the 
AAT and the SemelPol polarimeter. Applying the Least Squares 
Deconvolution method, they detected a magnetic field of around 50 
G in the Herbig Ae star HD 104237, but the y obtained no detec- 
tion i n another Herbig Ae/Be star, HD 1 00546. iHubrig et al.l J20041 
120071) claim detections in HD 139614 and HD 144432. Finally re- 
cent search for m agnetic field in the field Herbig Ae/Be stars by 
IWade et al] d2007l) has been carried out using the FORS1 spec- 
tropolarimeter at the ESO VLT. They identified two possible mag- 
netic stars (HD 101412 and BF Ori), of which one has been con- 
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firmed (HD 101412), with higher resolution data (Wade et al. 2007 
in preparation). 

Recently, the new generation instrument ESPaDOnS was in- 
stalled at the Cassegrain focus of the CFHT (Donati et al., in prepa- 
ration). Thanks to the high efficiency and qualities of this high res- 
olution spectropolarimeter, magnetic fields in several Herbig stars 
have been discovered. During the technical run of ESPaDOnS a 
magnetic field was discovered in the Herbig Ae/Be star HD 200775 
(Donati et al., in preparation). Then, during scientific nights of the 
first semester of ESPaDOnS, three other stars were detected as 
magnetic: HP 72 106, V380 Ori and HD 190073 dWade et al.l20ol ; 
ICatala et al.ll2007h . These recent discoveries bring new arguments 
to support the fossil field theory. However further exploration is 
necessary. 

According to the fossil field hypothesis, the structure and the 
intensity of the magnetic field of these PMS stars should be similar 
to those of the main sequence magnetic stars. In order to verify 
this, we have observed HD 200775, with ESPaDOnS, during many 
nights. Using the temporal variations of the Stokes V profile we 
were able to model the structure and the geometry of the magnetic 
field. The next section presents the observations and data reduction 
procedures. In Sect. 3, 4 and 5 we present the study of the intensity 
spectrum and in Section 6 the study of the polarization spectrum. 
Conclusions are given in Sect. 7. 



2 OBSERVATIONS AND DATA REDUCTION 

Our data were obtained using the high resolution spectropolarime- 
ter ESPaDOnS installed on the 3.6 m Canada-France-Hawaii Tele- 
scope (Donati et al., in preparation) during many scientific runs. 
Table[T]presents the log of the observations. One of our spectra has 
also been obtained using the instrument Narval, which is a copy of 
ESPaDOnS, installed on the 2 m Bernard Lyot Telescope (TBL) at 
the Pic du Midi observatory in France. 

We used the instruments ESPaDOnS and Narval in polarimet- 
ric mode, and we obtained spectra with resolving power of 65000. 
Each exposure was divided in 4 sub-exposures of equal time in or- 
der to compute the optimal extraction of the polarisation spectra 
(Donati et al. 1997, Donati et al., in preparation). We recorded only 
circular polarisation, as the Zeeman signature expected in linear po- 
larisation is about one order of magnitude lower than circular po- 
larisation. The data were reduced using the "Libre ESpRIT" pack- 
age especially developed for ESPaDOnS and Narval, and installed 
at the CFHT and at the TBL (Donati et al. 1997, Donati et al, in 
preparation). After reduction, we obtained the intensity Stokes / 
and the circular polarisation Stokes V spectra of the star observed, 
both normalised to the continuum intensity of HD 200775. A null 
spectrum (N) is also computed in order to diagnose spurious po- 
larisation signatures, and to help to verify that the signatures in the 
Stokes V spectrum are real jPonati et al.ll 19971) . 

The data of May 2005 were affected by a 1.3 mag loss com- 
pared to the data obtained more recently with ESPaDOnS. This 
problem, which was due to damage to the external jacket of the op- 
tical fibres, was fixed prior to the July run. In order to compensate 
for this damage and to obtain a satisfactory S/N ratio, we exposed 
longer (see Table [TJ. 

We then applied the Least Squares Deconvolution procedure 
to all spectra jPonati et alj 19971) . This method assumes that all se- 
lected lines of the intensity spectrum have a profile of similar shape. 
Hence, this supposes that all lines are broadened in the same way. 
We can therefore consider that the observed spectrum is a convolu- 



tion between a profile (which is the same for all lines) and a mask 
including all choosen lines of the spectrum. We therefore apply a 
deconvolution to the observed spectrum using the mask, in order 
to obtain the average photospheric profiles of Stokes / and V. In 
this procedure, each line is weighted by its signal to noise ratio, 
its depth in the unbroadened model and its Lande fa ctor. The mask 
was first computed using Kurucz ATLAS 9 models jKurucz|[T993rj 
with T e g = 19000 K and logg = 3.5 suitable for the star (Table 
We excluded from this mask hydrogen Balmer lines, strong res- 
onance lines and lines whose Lande factor is unknown. Then we 
cleaned the mask by selecting the lines with a depth less than 0.4, 
in order to eliminate the lines contaminated by strong emission. 
We have also modified the line depths in order to take into account 
the relative depth of the lines of the observed spectrum. The final 
mask contains only 37 lines, but is sufficent to compute the LSP 
Stokes / and V profiles with a good S/N ratio (Table [T). To ex- 
plore the sensitivity of the LSD profiles to the detailed line list, 
line masks were constructed using weak lines including He lines, 
and using weak lines excluding He lines. The various LSD profiles 
were analysed in the manner described in Section 3. No significant 
differences were found. The null N profile has been computed in 
the same way, from the null spectrum, and will be used, in Sect. 6, 
to check that the signature in the Stokes V profile is real. 

The LSD average line profiles were computed on a velocity 
grid with 3.6 km.s~' sampling. The resulting relative noise in the 
LSD Stokes V profiles is given in the 5th column of TableQ] 



3 BINARITY 

In the litterature, many autho rs refer to the binarity of HD 200775. 
ICorporon & Lagrange! dl999h detected radial velocity variations of 
photospheric lines of the spectrum of HD 20 0775, but they were no t 
able to determine the orbital period. iMiroshnichenko et alj d 19981) 
analysed spectroscopic data of this star obtained by many authors 
over 20 years. They noted a cyclic variation of the equivalent width 
of Ho- with a 1345 day period, and they sug gested that a compan- 
ion star could be the trigger of the Ha activity. iPogodin et al. I d2004h 
discussed again this theory using new data and observational evi- 
dence of radial velocity variations of some photospheric lines. They 
plotted two radial velocity curves, obtained on the one hand from 
the Ho- wings, and on the other hand from photospheric lines. They 
fitted them with a synthetic radial velocity curve in order to ob- 
tain the ephemeris and the or bital parameters of the double system. 
Finally iMonnier et alj d2006t) were able to separate both stars of 
HD 200775 using interferometry, and determine orbital parameters 
for the system. 

Figure [T] shows the LSD average Stokes / profiles (thin full 
line) of all our observations, obtained using the cleaned mask as 
explained in Section 2. We can see the presence of a second com- 
ponent in each profile : on the blue for the Sept. 2004 profile and on 
the red for the others. Furthermore, Fig. [2] shows some individual 
absorption lines of the intensity spectrum, observed in Sept. 2004 
and May 2005, where the secondary component is visible. Between 
Sept. 2004 and May 2005 the primary component of the profile and 
the photospheric lines shifted towards the blue, while the secondary 
component moved in the other direction. We therefore assign this 
second component of the average profile to a companion star. HD 
200775 is therefore an SB2 system whose orbital period is likely 
greater than one year. In the following we call HD 200775A (with 
the shaper, deeper lines) the primary component and HD 200775B 
(with the broader, shallower lines) the secondary component. 
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Figure 1. LSD Stokes / profiles of all spectra (thin full line) superimposed with their fit (thick dashed-line) described in Section 3. 
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Figure 2. Photospheric lines of Sill 4128A, SII 4153A, SII 5032A, and 
SII 5454A, oberved on September 23rd 2004 (redshifted) and on August 
25th 2005 (blueshifted). The fit (thick dashed line) with a convolution of 
a gaussian and a rotation profile is superimposed to the observed profiles 
(thin full line). 



In order to study the variations of the radial velocities of our 
spectra, we performed a simultaneous Least-Squares fit to the 30 
LSD / profilesQ Each profile is fitted with the sum of two func- 



tions, each function modeling the line profile of one components. 
Each one of these two functions is the convolution of a rotation 
function (for which the projected rotational velocity v sin i is a free 
parameter in the fit) and a gaussian whose width is fixed and com- 
puted fro m the spect ral resolution and the inferred macroturbulent 
velocity jGravll 19921) . We adopted an isotropic macroturbulent ve- 
locity of 15 km.s -1 in order to fit the wings of the profiles of in- 
dividual lines in the spectrum (both strong and weak lines), and of 
the LSD / profiles. 

The free parameters of the fitting procedure are the centroids, 
depths, and projected rotational velocities (v sin z) of both compo- 
nents. The centroids of both functions can vary from one profile to 
another., whereas the depths and vsinz' of both components can- 
not. This fitting procedure therefore assumes that the depths and 
v sin i of both components do not vary with time, which we con- 
firm (within the error bars) by fitting each profile separately. Figure 
[3] shows an example of such a fit for the Sept. 2004 and 22 May 
2005 profiles, as well as the profiles of the isolated components 
calculated by subtracting the synthetic profile of the other compo- 
nent from the observed LSD Stokes / profile. This automatic fitting 
procedure enables us to measure the v sin i and radial velocities of 
both components. We obtain projected rotational velocities of 26±2 
km.s" 1 and 59 ± 5 km.s -1 for the primary and the secondary com- 
ponents, respectively. Assuming an orbital period greater than one 
year (consistent with our data) and considering a projected rota- 
tional velocity of the primary of 26 km.s -1 , we conclude that the 
system is not synchronised. 

The inferred macroturbulent velocity is not determined very 
precisely: acceptable values range from about 10-30 km.s -1 . How- 
ever, the fitting procedure is largely insensitive to these uncertain- 
ties - changing the macroturbulence within the uncertainties results 



Although a potentially more effective approach would be to perform spec- 
tral disentangling, our data set was not suitable for such an analysis due to 



the sparse coverage of the orbital period and the small velocity separation 
of the lines of the components (O. Kochukhov, private communication). 
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Figure 3. Upper panels: LSD Stokes / profiles from spectra obtained on 
September 23, 2004 and September 21, 2005 (full line). The fitted line pro- 
files are the superimposed dashed lines. Lower panels: Individual intensity 
profiles of the primary (narrow red profile) and secondary (large green pro- 
file) components of the system. 



Figure 4. Upper panel : Equivalent width of Ha in function of time. Ha ap- 
pears in two different orders in the spectrum. Therefore the values measured 
on both orders are plotted. Lower panel : Radial velocities of the primary 
(red filled circle) and the secondary (green filled diamond) components of 
HD 200775. The full lines are the fitting curves. The arrows indicate the 
periastron passage. 



in only very small changes to the inferred v sin is, and to no signif- 
icant change in the radial velocities. 

We observe a strong variation of the radial velocity of the pho- 
tospheric lines of the primary star between Sept. 2004 and May 
2005 from — h7 km.s~' to 20 km.s~'. Then it is slightly vary- 
ing, and since August 2005 it has been increasing slowly. In Fig. [4] 
(lower panel), we plot the radial velocity of both components as a 
function of time, fitted with the radial velocity curves of an eccen- 
tric binary system, using a^ 2 minimization method. We obtain the 
following parameters for the system: orbital period P = 1412 ± 54 
d, periastron epoch To = 2448991 ± 152 d, systemic radial veloc- 
ity y = -7.9 ± 0.9 km.s _1 , eccentricity e = 0.32 ± 0.06, periastron 
longitude a> = 216 ± 12°, and radial velocity amplitude of both 
components Kp = 20.9 ± 2.5 km.s~' and K$ = 17.0 ± 2.5 km.s" 1 , 
leading to the mass ratio of the system q = = || = 0.81 ± 0.22. 

The orbital parame t ers th at we found are similar to those ob- 
tained by Pogo din et alJ d2004l) from the bisector velocities of Ho-. 
In order to compare our work to theirs, we also measured the bisec- 
tor velocities of the Ho- line, as well as the other emission lines of 
the intensity spectrum. We observe that all of them closely follow 
the radial velocity of the secondary star (see Fig. [5}- We therefore 
conclude that the region from which these emission lines originate 
is linked to the secondary component of the system. These results 
are in contradiction with those of Pogodin et al. who found that 
the bisector velocity follows the radial velocity of the primary star. 
Whether this contradiction is related to a real change of behavior of 
the binary system and its environment, or to problems of measure- 
ment of bisector and/or photospheric radial velocities in the data of 
Pogodin et al., remains to be determined. 

Assuming that the variations of the bisector velocities of Ha 
originate from the orbital movement of the secondary, we included 
them into the fit of the radial velocities. We found no significant 
differences of the inferred orbital parameters compared to those ob- 
tained from the fit of solely the radial velocities. 
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Figure 5. Radial velocities of the primary (red filled circle) and the sec- 
ondary components (green filled diamond) of HD 200775 superimposed 
with the bisector velocities of the Ha line (black open square), measured 
at a level of 1,5 — 2.0F C . Note that the bisector velocities trace the radial 
velocities of the secondary, and clearly do not follow those of the primary. 



Using our determination of the orbital parameters we can de- 
rive the semi-major axis of the orbit, and the masses of the pri- 
mary and the secondary, as functions of the inclination of the or- 
bit: a sin! orb , M P sin 3 i orb , and M s sin 3 i orb (see Table[2j. Using our 
determination of M P and M s (Table [3} from stellar evolutionary 
models, we can therefore estimate the orbital inclination of the sys- 
tem. We find sin i orb (P) = 0.70 ± 0. 16 from M P and M P sin 3 ! OTb , and 
sin / orb (S) = 0.78 ± 0. 17 from M s and M s sin 3 i orb . We determine the 
final value of sin ; orb by taking a weighted avarage of sin i or b(P) and 
sini orb (S). We obtain sini orb = 0.74 + 0.17, leading to i OTb = 48*!Z °. 
Then we use this value to derive the semi-major axis of the orbit 
: a = 6.7 ± 1.9 AU. Using the Hipparcos parallax of HD 200775 
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Table 2. Orbital parameters of the system found in this work, compared to previous work (Pogodin et al. 2004; Monnier et al. 2006). 



Parameter Radial Velocity Bisector velocity Interferometry 

This work (Pogodin et al. 2004) (Monnier et al. 2006) 



Period (days) 
To (HDJ) 

y(km/s) 
e 

o>n 

K P (km/s) 
(km/s) 



Results of RV fit 
1412 ±54 
2448991 ± 152 
-7.9 ± 0.9 
0.32 ± 0.06 
216 ± 12 
20.9 ± 2.5 
17.0 ±2.5 



1341 ±23 
2449149 ± 87 
5.0 ±0.6 
0.29 ± 0.07 

203 ± 22 

11.2 ±0.7 



1377 ±25 
2449152 ± 90 

0.30 ± 0.06 
224 ± 16 



a sin i m -b (AU) 
cip sin i orb (AU) 
«s sin i orb (AU) 



M s sin 3 r orb (M Q ) 

'orb (°) 
a (AU) 
a (mas) 



Derived parameters 
4.9 ±0.8 
2.73 ± 0.46 
2.21 ±0.38 
3.6 ± 1.6 

4.4 ± 1.9 

48+17 



6.7: 
16: 



1.9 



1.38 ±0.11 



65 ±8 
15.14 ±0.70 



(n = 2.33 + 0.62 mas), this corresponds to a projected separation of 
the components a = 16 ± 9 mas. 

The values of the orbital parameters that we derived are sum- 
marised in Tableland compared to the orbital elements of Pogodin 
et al. (2004), as well a s to the orbital param eters determined from 
interferometric data by Monnier et al. ( 2006). Most of our measure- 
ments of the orbital parameters agree well with both works. Our 
determination of A"s is larger than that of Pogodin et al. which is 
likely due to our lack of data during the minimum separation of 
both stars. We need to continue to observe the system in order to 
confirm that value. 



4 FUNDAMENTAL PARAMETERS OF THE PRIMARY 
AND THE SECONDARY 

We used the fundam ental parameters of HD 200775 of 
Herna ndez et alj J2004l) . These authors studied the spectra of 75 
stars, most of which are classified as Herbig Ae/Be stars. They 
determined the spectral type of each star as well as the luminos- 
ity with particular attention to the total-to-selective extinction (R v ) 
used. In the case of classical main sequence stars the interstellar 
extinction law gives R v = 3.1. However, the circumstellar mat- 
ter of Herbig stars may be dominated by dust and can lead to a 
higher reddening of the star . Using UBVR photome tric data of 
iHerbst & Shevchenkd ( 1 19991) . iHernandez etal] fc004l) found that 
Rv = 5 fits better most stars of their sample than the standard value 
R v = 3.1. The knowledge of R v is fundamental to determine the 
luminosity: using R v = 3.1, they found log(L/L ) = 3.73 ± 0.27, 
while with R v = 5, they found log(L/L Q ) = 4.17 ± 0.27. For the 
purposes of this analysis, we adopt the most likely value of R v = 5, 
as HD 200775, illuminating the reflect ion nebulae NGC 7 023, is 
still largely surrounded by dust and gas jFuente et al.ll 1998h . which 

certainly have an impact on the extinction law for this star. 

The observed luminosity, derived by Herna ndez et alj d2004l) . 
corresponds the total of the luminosities of both stars. A priori, we 
do not know the luminosity ratio of the components. However, we 
can see in the spectrum that the lines of the secondary component 
are less deep than those of the primary. We suppose that this is 



mainly due to its higher rotational velocity (vsin/ = 59 km.s" 1 , 
Section 3) and slightly due to a fainter luminosity. Furthermore, all 
lines observed in the spectrum of the secondary are also observed in 
the spectrum of the primary star. For all these reasons, we suspect 
that the temperature of the secondary star is similar to that of the 
primary component. 

Hernandez et al. (2004) determined the spectral type of 
HD 200775 mainly by comparing the strength of atomic absorption 
lines to those of standard stars and found log(r c ff) = 4.27. We com- 
pared our spectra with a synthetic spectrum of a double star whose 
components have the same effective temperature and surface grav- 
ity : T e ff = 19000 K and logg = 3.5. First, we calculated separate 
spectra of both components using TLUSTY non - LTE atmosphere 
models and the SYNSPEC code dHubenvl[l988l: iHubenv & Land 
Il992lll995l) . Then we computed the spectrum of the double star, 
with the BINMAG1 program (O. Kochukhov, private communica- 
tion), using a macroturbulent velocity of 15 km/s, our measured 
v sin i, and radial velocities, and our estimation of the luminosity 
ratio L s /Lp = 0.67 (see below). We conclude that the temperatures 
of the synthetic spectrum are in good agreement with the observed 
spectrum of HD 200775. Figure[6]shows a portion of the observed 
spectrum superimposed on the synthetic one. 

In addition, we also compared our spectra with a synthetic 
spectrum of a single star of r cff = 19000 K and logg = 3.5. We 
found that the synthetic spectrum of a double star fits better the 
depth and the width of the lines than the spectrum of a single star. 
That adds additional support to our hypothesis that both compo- 
nents of the system have similar temperature. In order to quantify 
the temperature range of the stars, we have calculated synthetic 
spectra by varying the effective temperature of both stars around 
19000K. We found that these spectra are able to fit our observa- 
tions within ±2000 K around 19000K. 

The luminosity of Hernand ez et alj J2004h is the sum of the 
luminosity of the primary (L P ) and the secondary (L s ) components 
of the system. We can estimate the luminosity ratio rt = Ls/Lp 
from the Stokes / profiles (normalised to the continuum of the bi- 
nary). Section 3 describes the fitting procedure of these profiles. 
The result of the fitting procedure is the function: 
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Figure 6. ESPaDOnS spectrum (full line) of HD 200775 between 438 and 450 nm supe rimposed with the synthetic spectrum of a double star (dashed-line) 
calculated using TLUSTY non-LTE atmosphere models and the SYNSPEC code faubenvfl 988; Hubeny & Lanz 1992l ll995h . 



Table 3. Fundamental parameters of HD 200775 A and HD 200775B. Column 7 gives the age of the system. 



Star 


Teg 
(K) 


log(L/L Q ) 


M/M Q 




(cm.s -2 ) 


age 
(Myr) 


v sin i 
(km.s- 1 ) 


HD 200775A 
HD 200775B 


18600 ±2000 
18600 ±2000 


3.95 ± 0.30 
3.77 ± 0.30 


10.7 ±2.5 
9.3 ±2.1 


10.4 ±4.9 
8.3 ± 3.9 


3.4 ± 1.2 
3.6 ± 1.2 


0.1 ±0.05 


26 ±2 
59 ±5 



/=1-./p-/s = 1- 



l+r L 



/p- 



I'L 



l+r L 



fs, 



(1) 



where / P and / s are the separate Stokes / profile shapes of the pri- 
mary and the secondary stars. The ratio J / s 'dv/ f fpdv is therefore 
equal to rt * W$/Wp, where Wp and Ws are the equivalent widths of 
the separate Stokes / profiles of both stars. These profiles have been 
computed using the same mask, and therefore the same lines. As- 
suming a common origin for both stars of the system and that they 
are not peculiar, they should have approximately the same chem- 
ical composition. The intrinsic equivalent widths should therefore 
be identical and the ratio j f$dv/ f fpdv be equal to r^. We com- 
puted this ratio for all our observations. The mean and the standard 
deviation of these 30 values gives a luminosity ratio L S /L P equal to 
0.67 ± 0.05. We therefore derive the luminosity of the stars using : 



log — 



log ■ 



log log(l + r L ) = 3.95 ±0.30 

log — - log(l + — ) = 3.77 ± 0.30 



(2) 
(3) 



where, L = L P + L s is the observed luminosity of HD 200775 
l lHernandez et al.ll2004[) . 

We used the derived luminosities and the temperature of HD 
200775 to place both stars in the HR diagram and to compare their 



position with evolutionary tracks of differe nt masses ca lculated 
with the CESAM stellar evolutionary code jMorel|[l997h . In this 
way we obtained the mass and the radius of both stars, as well as 
the age of the system, assuming a common origin for both stars (the 
age of the system is the intersection of both age ranges determined 
for each star separatly) (Fig.UJl. The value of the fundamental pa- 
rameters are summarized in Table [3] The mass ratio is therefore 
q = 1.1 ± 0.5, which is consistent with the value determined from 
the orbit analysis (q = 0.81 ± 0.22, Sec. 3). Note that in both cases 
the error bars on the mass ratio are very large. This is partly due to 
our poor coverage of the orbital phase : our data cover only the half 
of the orbital period. In particular we do not have data where the 
maximum radial velocities are predicted. Therefore we need to get 
more data during the next years to improve it. The large error bars 
of the orbital parameters are also due to the non-unique solution 
of the fit of the Stokes / profiles, leading to very large error bars 
on the radial velocities and the luminosity ratio. To better constrain 
the mass of the system, photometric observations of both separate 
components would be necessary. It would provide us the luminos- 
ity of each star and therefore we would no longer be dependent on 
our estimation of the luminosity ratio. 

Furthermore, we note that using the orbit analysis we find that 
the secondary is more massive than the primary, whereas in this 
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Figure 7. Evolutionnary tracks (full lines) of different mass : from 7 M Q 
to 16 M Q , plotted in a HR diagram. The dashed lines are the 0.01, 0.025, 
0.04 Myr, 0.065 Myr, 0.09 Myr, 0.115 Myr, 0.14 Myr, 0.17 and 0.2 Myr 
isochrones. Crosses represent the error bars in temperature and luminosity 
of HD 200775A (red square) and HD 200775B (green triangle). 
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Figure 8. Stokes / (bottom) and V (top) mean profiles of HD 200775 ob- 
served on Sept. 23rd 2004 (full line) and May 23rd 2005 (dashed line). The 
null profile N is plotted in the middle, and the 3cr error bars of V and N are 
plotted next to the respective profiles. 



analysis we find that it is less massive than the primary. This in- 
consistency could result from an inaccurate estimation of the lu- 
minosity of both components, determined by assuming a similar 
reddening for both stars. This uncertainty, which we are unable to 
take into account in our error bars, could lead to an underestima- 
tion of the luminosity of the secondary relative to the primary and 
therefore to an overestimation of the mass of the secondary relative 
to the primary. This point is discussed in detail in Sec. 7. 



5 PROPERTIES OF THE INTENSITY SPECTRUM 

5.1 Abundances anomalies 

Magnetic intermediate mass stars on the main sequence show 
strong photospheric abundance peculiarities. As some Herbig 
Ae/Be stars must be the evolutionary progenitors of these chemi- 
cally peculiar Ap/Bp stars, we searched for abundances anomalies 
in our spectra, as well as time profile variations due to abundance 
spots. We chose in the same spectral region lines of different chem- 
ical species whose depths predicted by the synthetic spectra are 
identical. Then in the observed spectra, we compared the depths 
of these lines. We used this method, insensitive to the veiling phe- 
nomenon, to try to find over- or under- abundant species. We found 
no systematic differences between observed and calculated equiva- 
lent widths of numerous chemical species. We conclude that there 
are no strong abundance anomalies, in the limit of ±0.40 dex in 
abundance, in the spectrum of HD 200775. In particular, we ob- 
serve no peculiarity in the helium lines nor any variability of the 
equivalent width of various species from one night to another and 
over many nights from Sept 2004 to April 2007. 

These results are different from the behavior of most of the 
hot magnetic B stars on the main sequence, which are frequently 
He-rich or He-weak, and which show abundances starspots varying 
with the rotational phase of the star. This suggests that at the young 
age of HD 200775, abundance anomalies have not had sufficient 
time to develop, or that they are limited by ongoing mixing due to 
accretion and mass loss. 

5.2 The emission lines 

We observe two kinds of emission lines in the spectrum of 
HD 200775 : broad lines and narrow lines. Amongst the broad 



emission lines we find mainly the Balmer lines and the OI 7772 A 
lines, and some Fe II and Si II with FWHM around 130 km.s~' . The 
quality of our data as well as the complex structure of the emission 
lines do not allow us to fit them with a gaussian function, and there- 
fore to measure the radial velocities of these lines. However the 
bisector velocities of these lines follow the radial velocity of the 
secondary (Sec. 3). We note a decrease of the equivalent widths of 
these emission lines from Sept 2004 to April 2007 (see Fig.|4]upper 
panel, for Ho-). We observe in Fig.[4]that the maximum equivalent 
width of Her occurs close to the periastron passage of the orbit (in- 
dicated by the arrows) that we derive i n Sect ion 3. This is consistent 
with the hypothesis o f |Pogodinetai] ( l2004l) who suggest that bina- 
rity could be the origin of these variations. 

In addition to these broad lines, we note a large number of nar- 
row emission lines, about a hundred over the whole spectrum, with 
a FWHM around 19 km.s -1 . Some of these are clearly visible in 
Fig. [6] We do not observe any shift in radial velocity of these lines 
during our 2- year observations . Using the numerous observations 
of t] Carinae dThackeravl [l953l 1 1 962l [l 9671 ; lHamann et all 1 19941) 
we were able to identify all of them and to conclude that they are 
nebular lines (lines coming from the nebulosity illuminated by HD 
200775). Appendix A gives the identification of all lines. 



6 MAGNETIC FIELD ANALYSIS 

Figure [8] shows the Stokes / and V LSD profiles obtained on Sept. 
2004 and May 23rd 2005. First we note that the null profile N (Sec. 
2) is totally flat indicating that the signature in the V profile is real. 
Secondly the V profile has shifted from the blue to the red, as the 
primary component of the / profile, between Sept. 2004 and May 
2005. This behaviour is observed throughout the entire data set. 
For this reason we attribute the observed magnetic field solely to 
the primary component. 

To investigate the possible presence of a magnetic field in the 
secondary, we calculated synthetic Stokes V profiles for a star sim- 
ilar to the secondary one, with a vsin/ of 59 km.s~', and consid- 
ering the same magnetic proprieties of the primary, described in 
Section 6.2. The maximum signal in circular polarisation predicted 
is around 10~ 3 , which is the same order of magnitude as the error 
bars of our observed V profiles. Therefore, if the secondary com- 
ponent of the system hosts a magnetic field of similar topology and 



Characterisation of the magnetic field of the Herbig Be star HD 200775 9 



similar intensity to that of the primary, it would be difficult to detect 
it with our data. On the other hand, if our polarisation spectra con- 
tain a signal of the secondary star corresponding to a similar field 
intensity, it is negligible with respect to the signal of the primary 
star. 



6.1 Phased longitudinal field variation 

Figure [TO] shows the variations of the 30 Stokes V profiles, nor- 
malised to the continuum intensity of the binary, observed from 
2004 to 2007. The longitudinal magnetic field B f , which is the pro- 
jected magnetic field onto the line of sight, can be o btained from 
both L SD S tokes I and V pro files, as explained by iDonati et al.l 
d 19971) and IWade et alj feOOd) . With this aim we compute the 
Stokes / and V profiles of the primary (Ip and Vp) as follows. Using 
Eq. (1) and the fit of the Stokes / profiles of the binary (Sec. 3), we 
subtract from the observed Stokes / profile the synthetic profile of 
the secondary component. Then we renormalised to the continuum 
of the primary: 



l-(l + ft)(l -I-f$l 



(4) 



As we assume that the observed magnetic field is only from the 
primary, we only need to renormalise the observed Stokes V profile 
by the continuum of the primary, in order to get Vp: 

V P = (1 + r L )V (5) 

Using these new profiles and Eq. (1) of lWadeetai] j2000), wc 
measured the longitudinal magnetic field B[ of our data. We fit the 
variations of B[ with a sinusoidal function of 4 parameters : 

• P, the rotation period of the star, 

• to, the reference Julian Date of the maximum of Be, 

• B, the semi-amplitude of the curve, 

• B , the shift of the sinusoidal curve with respect to B f = 0, 

consistent with a dipole centered inside the star. The best fit su- 
perimposed on the data in Fig. [9] gives the following parameters : 
P = 4.328 + 0.003 d, t a = 2453515.8 + 0.3 d, B = -309 ± 115 G 
and S = - 1 39 + 95 G, with a reduced x 1 = 1 • 1 • We used that value 
as a first estimation of the rotation period in the fitting procedure of 
the Stokes V profiles. 

This analysis reveals the basic properties of the magnetic field 
of HD 200775A: that it is organised on large scales; that it has an 
important global dipole component; and that the dipole geometry 
is stable on timescales of several years. 



6.2 Fitting of Stokes V profiles 

To fit the Stokes V profiles directly , we use th e dece ntred dipole 
oblique rota tor model describe d by lLandstree t (1970). Using the 
relations of lLandstreetl Jl970t) . giving the intensity of the mag- 
netic field inside a star, we calculate the longitudinal magnetic field 
be (6, tp), at each point {6,ip) of the surface of the star of velocity v, 
in classical spherical coordinates in the observer's frame. 

We assume a gaussian local intensity profile of width <x 
and depth d. The width is calculated using the resolving power 
of the instrument and the macroturbulent velocity determined in 
Section 3. The depth is determined by fitting the Stokes / pro- 
files of the primary component, determined in Section 6.1. We 
then calculated the local Stokes V profile at each point on the 
surface of the star using the weak magnetic fie ld relation of 
Landi deg lTnnocenti & Landi degl'Innocentil dl973h : 




Figure 9. Longitudinal magnetic field plotted in function of the rotation 
phase. The dashed line is the best sinusoidal fit. 



Table 4. Ranges and minimum bins of parameters explored in the fit of the 
Stokes V profiles 



Parameters 


Min 


Max 


bin 


P (days) 


4.0 


5.0 


0.0001 


T (HJD) 


2453514 


2453519 


0.02 


if) 





90 


1 


PC) 





180 


1 


Bd(G) 





2000 


10 


rfdip («,) 


-0.4 


0.4 


0.01 



V(y,6,tp) = -CgAocbc(6,ip) 



dv' 



(6) 



where C = 4.67 x 10~ 13 A~ 1 G"' 1 , g and A (A) are the mean Lande 
factor and wavelength of the lines used in the mask (Section 2), 
c is the speed of light, / is the local intensity profile, and v is the 
velocity. Then we integrated over the visible stellar surface usin g 
the limb darkening law with a parameter equal to 0.4 faarette OOO). 
We obtain the synthetic Stokes V(y) profile, that we normalised to 
the intensity continuum, to compare to the observed V profiles. This 
model depends on five parameters : 

• P the rotation period 

• t , the reference Julian Date of the maximum of the surface 
magnetic intensity, used with P to compute the rotation phase 

• i, the inclination of the stellar rotation axis to the observers 
line-of-sight, 

• B the magnetic obliquity angle, 

• B A the dipole magnetic intensity, 

• rfdip the displacement of the dipole from the center of the star, 
along the magnetic axis, in stellar radii (R,). 

We calculated a grid of V profiles for each date of observation 
(see Table|4]for details on the grid), varying the six parameters (and 
assuming for the initial value of P the solution obtained from mod- 
elling the longitudinal field variation). Then we applied &x' min- 
imisation to find the best model which matches simultaneously our 
30 observed profiles. The best model that we found, with^ 2 = 1.0, 
corresponds to P = 4.3281 ± 0.0010 d, i = 60 ± 1 1°, B = 125 ± 8°, 
S d = 1000 ± 150 G, and d dip = 0.05 ± 0.04 R„ where the error 
bars correspond to 3tr confidence. Figure [TO] shows the synthetic 
Stokes V profiles superimposed on the observed ones. We see that 
this model acceptably reproduces most of the observed V profiles. 
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Figure 10. LSD V profiles (noisy black line) of the 30 spectra superimposed to the best oblique rotator model (smooth red line). The numbers close to the 

profiles are the rotation phase, and the little bars on the right of the profiles are the mean error bars in V. The profiles are sorted by increasing date, as in Fig. 1. 



The small value of the dipole decentring parameter is nearly con- 
sistent with zero. The available data are therefore consistent with a 
dipolar magnetic field at the centre of the star 

In order to compare this result to the fit of the longitudinal field 
variations, we estimate the intensity and the obliquity angle /? 
of the magnetic dipole component using the fit of the longitudinal 
field variations. Using the well-kn own equations relating i, f$ and 
S(i to the longitudinal field extrema dBorra & Landstreetll980h . and 
adopting i = 60 ± 1 1° , we findfl = 128 ± 13° and B d = 1500!™ G. 
These values agree within the error bars with those determined with 
the fit of the Stokes V parameters, which confirms the topology 
of the magnetic field found using this fit. Table |5] summarises the 
parameters of the adopted dipole model. 

Because the Stokes / and V profiles are influenced in similar 
ways by uncertainties in the procedure that we have used to recover 
the individual LSD profiles of the primary, the magnetic modeling 
is only weakly sensitive to the details of the procedures described 
in Sect. 4 and Sect. 6. From a variety of experiments we find that 
neither the assumed luminosity ratio, nor the details of the fitting 
procedure yielding the v sin is, profile depths and radial velocities, 
have any significant direct influence on the derived magnetic mod- 
els. 

Finally, we point out that the magnetic modelling provides a 
strong constraint on the inclination of the rotation axis i. In partic- 
ular, for smaller values of i, we are unable to reproduce the Stokes 
V profile shapes. 



7 DISCUSSION 

7.1 Consequences of the magnetic model on the fundamental 
parameters of the star 

Using our determinations of the v sin i (from the LSD / profiles) 
and the period and inclination (from the magnetic model), we can 
estimate the radius of HD 200775A. We find R = 2.6 ± 0.5B o , 



Table 5. Magnetic dipole model of HD 200775A 



P(days) 4.3281 ±0.0010 

7/ (HJD) 2453515.8 ±0.8 

i(°) 60 ±11 

PC) 125 ±8 

fi d (G) 1000 ±150 

ddip (R«) 0.05 ± 0.04 



which is significantly smaller than the radius found from stellar 
evolutionary models (Table 0. In particular, this value suggest a 
star close to the main sequence, contrary to the stellar models. The 
model value is determined from the stellar effective temperature 
(which we consider to be well-determined) and the inferred lumi- 
nosity of the primary star. The primary's luminosity is itself deter- 
mined based on the ratio of equivalent widths of the LSD profiles, 
and from the dereddened photom etry of the combined system sup- 
ported bv lHernandez et alj (2004), which assumed equal reddening 
for each of the two components. 

The inferred luminosities could be affected by two possibly 
incorrect assumptions inherent in this analysis. First, based on our 
derivation of nearly identical effective temperatures for the two 
components, we have assumed that the intrinsic equivalent widths 
of their LSD profiles are equal. We then used the observed equiva- 
lent widths of the two components, determined from the LSD pro- 
files, to infer their relative luminosities. However, if the tempera- 
tures of the primary and secondary differ by the maximum allowed 
by the error bars (about 4000 K), the combined effects of differ- 
ences in excitation and ioniosation, and the T* s dependence of the 
flux contribution, could combine to generate important relative dif- 
ferences (either larger or smaller depending on the effective temper- 
atures) in the contributions of the spectra of the two components to 
the observed LSD profiles. 

The second assumption that could influence the inferred lu- 
minosities is that the two components are affected by similar red- 
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dening. In the i r inter ferometric study of the HD 200775 system, 
Mon nier et alj |2006) resolved the circumstellar disc of the sec- 
ondary and found that the secondary was brighter than the primary 
in H band by a factor of 6.5. If this IR flux difference results from 
the presence of CS material around the secondary that is not present 
around the primary (which would be consistent with our spectro- 
scopic observations, indicating that the Ho- emission is associated 
with the secondary), it would imply that the secondary is affected 
by additional reddening that has not been taken into account in 
our analysis. The increased extinction of the secondary would re- 
sult in an underestimated luminosity of this star relative to the pri- 
mary. Consequently, the primary's luminosity would be overesti- 
mated, leading to a systematic overestimation of its inferred radius. 
It would also lead to an overestimation of the mass of the primary 
relative to the secondary and therefore to a mass ratio q = M P /M S 
greater than one (as observed), which is inconsistent with the orbit 
determination of q = 0.81 ± 0.22. 

Neither of these potentially important sources of systematic 
error can be investigated in detail with the available data, and ul- 
timately we are led to conclude that the primary's radius could be 
anywhere from the ZAMS radius to the radius reported in Table 3. 
Likewise the masses and ages of the stars derived in Sec. 4 from 
the luminosity and temperature must be considered with caution. 
Individual photometric observations of each of the stars is required 
to get independent measurements of the luminosities of both stars, 
and therefore accurate values of their masses, radii and ages. 

7.2 Consequences of the magnetic model on the origin of the 
system 

Assuming i = 60° implies that the primary would be close to the 
main sequence, and therefore would have a lower luminosity than 
the secondary. According to the luminosity ratio in H band of Mon- 
nier et al. (2006), the secondary is redder that the primary and there- 
fore seems to be surrounded by denser circumstellar matter than the 
secondary, which is consistent with our spectroscopic data. This 
would indicate that the secondary is in a younger evolutionary state 
than the primary. Furthermore the eccentricity of the orbit found 
in Sec. 3 is very large (Table 2). These properties may well sug- 
gest that the two components are at two very different evolutionary 
stages, which could imply that the system formed by capture in- 
stead of forming initially as a simple close binary. 



8 CONCLUSIONS 

In the framework of the understanding of the origin of the magnetic 
field of the intermediate and high mass stars, we have begun to 
acquire spectropolarimetric observations of their potential progen- 
itors, the Herbig Ae/Be stars, including HD 200775. A magnetic 
field has been discovered in the latter and in order to determine its 
topology, we monitored it during many successive nights and over 
more than 2 years, using the high resolution spectropolarimeters 
ESPaDOnS and Narval. 

First we inspected the intensity spectrum and found that this 
star is a double-lined spectroscopic binary, whose secondary seems 
similar, in temperature and luminosity, to the primary. We mea- 
sured accurate values of the projected rotational velocity of both 
stars: 26 ± 2 km.s -1 and 59 ± 5 km.s -1 . We fitted the radial veloc- 
ity curves of both components and found an ephemeris similar to 
that of Pogodin et al. (2004). The amplitude of the radial velocity 
curves leads to a poorly constrained mass ratio. Observations of the 



separate components, in order to better determine the luminosity 
ratio of the system, are required to determine accurately the lumi- 
nosities and hence the masses of the components. According to our 
determination of the orbital period and of the masses of both stars, 
we estimate the separation between both components of the system 
aroun d 16 mas (the distance of the system is estimated around 430 
pc bv lvan den Ancker et al.1 dl997h ). Only interferometric observa- 
tions of the system would be able to separate the two components 
and to estimate the luminosity ratio. 

We have also shown that the broad emission observed in the 
intensity spectrum is linked to the secondary component of the sys- 
tem. We found no evidence of abundance anomalies in the spec- 
trum. Finally, many narrow emission lines have been observed over 
the whole spectrum, which we identify as nebular lines. 

The position and velocity variations of the Stokes V profile 
shows that the detected Zeeman signature corresponds solely to the 
primary component. However, if the secondary hosts a magnetic 
field of similar topology and intensity to the primary, the signal in 
polarisation would be negligible with respect to the signal of the 
primary. On the other hand, this implies that our observations are 
relatively insensitive to the presence of a field in the secondary. 

We modeled the temporal variations of the Stokes V profiles 
in two different ways. First we measured the longitudinal magnetic 
field of the star and we fitted a sinusoidal curve, as predicted by 
the dipole oblique rotator model. Then we used the period found as 
a first estimation of the rotation period of the star in the Stokes V 
profiles fitting procedure. We considered an oblique rotator model 
and we applied x 2 minimisation to match simultaneously the 30 
Stokes V profiles that we observed over more than one year. We 
find that this simple decentered dipole model is sufficient to fit most 
of the profiles. In this way we provide the magnetic field topology: 
a dipole of intensity ~ 1000 G, displaced by 0.05 ± 0.04 R, from 
the center of the star towards the positive magnetic pole, whose the 
rotation axis is inclined by 60 + 11° with respect to the observer 
line-of-sight, and magnetic axis is inclined by 125 ± 8° with respect 
to the rotation axis. HD 200775 A rotates with a period of 4.3281 ± 
0.0010 d. 

We therefore show that the magnetic field of this star is ap- 
proximately dipolar, strongly inclined with respect to the rota- 
tion axis, with a polar intensity of 1000 G, and stable over more 
than two years. These characteristics ar e similar to the topology of 
the magnetic fields of the Ap/Bp stars teorra & Landstreetj ll980l : 
lBohlenderetal.ll 19871) . 

We have concluded that significant uncertainties exist related 
to the luminosity, and therefore the mass, radius and age of the mag- 
netic primary star. To clarify these uncertainties and to move for- 
ward, additional photometric observations of the individual com- 
ponents are required. In addition, a more complete spectroscopic 
coverage of the orbital cycle is needed, with the hope of ultimately 
performing spectral disentangling to yield the individual spectra of 
the primary and secondary. 
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APPENDIX A: NEBULAR LINES 



Characterisation of the magnetic field of the Herbig Be star HD 200775 



Table Al: Identification of the nebular lines in the spectrum of HD 200775. Col- 
umn 1 gives the wavelength observed in our spectra, column 2 gives the equiv- 
alent width of the line, column 3 gives the ion, column 4 gives the multiplet 
number, colu mn 5 give s the re st wavelength of the l ine, and column 6 give s the 
reference. I: IThackeravl Jl953h . II: IThackeravl dl962h . Ill: IThackeravl dl967t) . IV: 
lHamann et alj dl994h . 



A (A) 


W A (mA) 


Ion 


Identification 


(A) 


Rf. 


4114.34 


8.5 


Fell 


23F 


4114.48 


Ill 


4177.07 


10.9 


Fell 


21F 


4177.21 


III 


4210.99 


PCI^ 


Fe II 


23F 


4211.10 


III 


4243.83 


42.6 


Fell 


21F 


4243.98 


III 


4244.68 


11.7 


Fell 


21F 


4244.81 


III 


4276.69 


40.7 


Fell 


21F 


4276.83 


I 


4287.25 


85.5 


Fell 


7F 


4287.40 


I 


4305.75 


9.3 


Fell 


21F 


4305.90 


1 


4319.47 


19.7 


Fell 


21F 


4319.62 


1 


4326.13 


5.4 


Nill 


3F 


4326.28 


1 


4346.72 


10.0 


Fell 


21F 


4246.85 


I 


4352.64 


19.5 


Fell 


21F 


4352.78 


I 


4358.23 


22.9 


Fell 


21F 


4358.37 


I 


4359.18 


70.6 


Fell 


7F 


4359.34 


I 


4372.28 


9.4 


Fell 


21F 


4372.43 


1 


4382.57 


5.2 


Fell 


6F 


4382.75 


1 


4413.62 


56.4 


Fell 


7F 


4413.78 


I 


4416.11 


51.1 


Fell 


6F 


4416.27 


I 


4451.95 


30.2 


FeH 


7F 


4452.95 


I 


4457.80 


29.7 


Fell 


6F 


4457.95 


I 


4474.75 


20.0 


Fell 


7F 


4474.91 


1 


4488.60 


12.1 


Fell 


6F 


4488.75 


I 


4492.49 


7.6 


Fell 


6F 


4492.64 


I 


4414.74 


7.4 


Fell 


6F 


4414.90 


I 


4528.23 


6.6 


Fell 


6F 


4528.39 


I 


4639.51 


12.9 


Fell 


4F 


4639.68 


I 


4664.29 


4.5 


Fell 


4F 


4664.45 


I 


4727.91 


24.9 


Fell 


4F 


4728.07 


I 


4773.56 


17.4 


Fell 


4F 


4774.74 


I 


4814.38 


54.2 


Fell 


20F 


4814.49 


I 


4874.33 


16.3 


Fen 


20F 


4874.49 


I 


4889.46 


33.6 


Fell 


4F 


4889.63 


I 


4898.44 


9.7 










4905.19 


29.5 


Fell 


20F 


4905.35 


1 


4947.23 


7.5 


Fell 


20F 


4947.38 




4950.58 


12.1 


FeH 


20F 


4950.74 




4973.22 


15.1 


Fell 


20F 


4973.39 




5005.37 


13.7 


Fell 


20F 


5005.52 




5020.08 


13.6 


Fell 


20F 


5020.24 




5040.81 


26.7 










5043.81 


6.0 


Fell 


20F 


5043.53 




5107.75 


7.6 


Fell 


18F 


5107.95 




5111.46 


17.7 


Fell 


19F 


5111.63 




5157.83 


17.9 


Fell 


18F 


5158.00 




5158.61 


53.5 


Fell 


19F 


5158.81 




5163.79 


20.8 


Fell 


35F 


5163.94 




5168.85^ 


9.8 


Fell 


42 


5169.03 




5181.78 - 


10.0 


Fell 


18F 


5181.97 




5199.00 


8.0 










5219.89 


17.7 


Fell 


19F 


5220.06 


III 



Inverse Pcygni profile, the equivalent width have therefore not been determined 
Blended with 5168.73 unidentified 
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Table Al: continued. 



X (A) 


W x (mA) 


Ion 


Identification 


(A) 


Rf. 


5261.46 


64.4 


Fell 


19F 


5261.61 


Ill 


5268.70 


15.6 


Fell 


18F 


5268.88 


III 


5273.18 


4.3 


Fell 


18F 


5273.38 


III 


5282.94 


5.9 










5296.67 


12.6 


Fell 


19F 


5296.84 


III 


5333.48 


39.7 


Fell 


19F 


5333.65 


III 


5347.50 


4.7 


Fell 


18F 


5347.67 


III 


5376.28 


32.8 


Fell 


19F 


5376.47 


III 


5412.49 


9.1 


Fell 


17F 


5412.64 


III 


5414.80 


4.6 










5432.99 


16.6 


Fell 


18F 


5433.15 


III 


5477.05 


9.4 


Fell 


34F 


5477.25 


III 


5527.17 


19.3 


Fell 


17F 


5527.33 


III 


5556.18 


2.1 


Fell 


18F 


5556.31 


III 


5580.65 


4.8 


Fell 


39F 


5580.82 


III 


5673.01 


8.1 


Fell 


F 


5673.22 


III 


5746.78 


22.0 


Fell 


34F 


5746.96 


III 


5835.25 


6.6 


Fell 


F 


5835.44 


III 


6364.94 


9.0 










6371.07 


24.1 


Sill 


2 


6371.36 


III 


6666.60 


13.4 


Mil 


2F 


6668.8 


III 


6729.64 


4.0 


Fell 


31F 


6729.85 


11 


6808.99 


6.0 


Fell 


31F 


6809.21 


II 


6813.41 


6.7 


Mil 


8F 


6813.73 


II 


7154.93 


51.8 


Fell 


14F 


7155.16 


IV 


7171.74 


18.2 


Fell 


14F 


7172.00 


IV 


7377.63 


73.5 


Mil 


2F 


7377.90 


IV 


7387.91 


15.9 


Fell 


14F 


7388.18 


IV 


7411.41 


31.4 


Mil 


2F 


7411.60 


IV 


7452.30 


24.5 


Fell 


14F 


7452.54 


IV 


7495.36 


11.3 


Fell 




7495.62 


IV 


7512.91 


25.6 


Fell 




7513.16 


IV 


7731.41 


8.3 


Fell 




7731.68 


IV 


7999.80 


129.4 


CrII 


IF 


7999.85 


IV 


8125.05 


104.0 


CrII 


IF 


8125.22 


IV 


8184.62 


2.7 










8216.04 


97.8 


NI 


2 


8216.28 


IV 


8308.22 


60.1 


CrII 


IF 


8308.39 


IV 


8616.65 


65.1 


Fell 


13F 


8616.96 


IV 


8891.62 


26.0 


Fell 


13F 


8891.88 


IV 


8229.43 


64.6 


CrII 


IF 


8229.81 


II 


8287.53 


56.3 










8694.94 


31.8 


SI 


6 


8694.70 


II 


8715.46 


11.1 


Fell 


42F 


8715.84 


II 


8768.89 


23.4 










8819.24 


17.2 










9122.64 


21.9 


Fell 




9122.91 


IV 


9186.94 


16.6 


Fell 




9187.15 


IV 


9203.80 


42.2 


Fell 




9204.05 


IV 


9217.98 


86.7 










9243.94 


177.7 
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